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Note on the Theory of Optical Images. 

By Geobge Steic. 



L. Mandelstam* obtains the result, that for a given aperture of a diaphragm 
only those linear structures will be similarly reproduced — through a lens-system 
— which satisfy the integral equation 

f(Z)=Z.f±Zf(x)K(x-Z)dx, 

in the case of an illuminated structure ; and the integral equation 

f\Z) = lS + -lf(<c)K*< < x-Z)dx, 

in the case of a self-luminous structure. Here f(x) denotes the distribution 

of amplitudes in the linear object, or — briefly — the structure itself. 

K(x — £)=K(s) is a function which depends only on the aperture of the 

J (c- s) 
diaphragm; when the aperture is circular, we have K(s) = -^ '- .f X is a 

s 

positive constant. 2 a is the length of the linear object ; the problem is more 
general, when the length of the structure is taken to be infinite. 
Mandelstam assumes for the solution of the above equations 

f(x) = cos(bx+q>), 

which is permissible on certain suppositions. For the reproduced images we 
have the relation 

where — satisfies, for illuminated and self-luminous structures of infinite 
length, the equations 

^- = 2fZK(s)'COsbs-ds 



and 
respectively. 



*6 



±=2fiK i (s)-cosbs-ds 



* Weber's Festschrift (1912), "ttber eine Anwend. d. Tntegralgleich.," p. 237. 

2 irr 
t J t (c-s) = Bessel's function; c== -—■ = const.; r = radius of the aperture; I = wave length of 



light; / = focus distance of the lens. 

13 
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Mandelstam obtains explicit results for — in the case of a square aperture 

of the diaphragm, both for illuminated and for self-luminous structures ; he ob- 
tains also an explicit result for a circular aperture, but only for illuminated 
structures. In the case of a circular aperture in the diaphragm and self- 
luminous structures of infinite length, he points out, with Lord Rayleigh,* that 
the solution of the problem depends on the evaluation of the integral 

f°°[ Jl{c s ' s) ~t- cos bs-ds, (I) 

which figures in the expression of r- for self-luminous structures and circular 
aperture of the diaphragm. We easily recognize that such functions as 
j- JxCcgn 2 . cosbs or r Vi(cs) T . sinbs> J^cs) = Bessel's function, 

are finite and continuous for all values of s. 

The purpose of this note is the evaluation of the former integral for the 

case b = 2c. 

We will start by taking into consideration the following integral : 

* (e) =C XX [^p 2 " dz ds, (II) 



where e = «- and i = V — 1, Ji(s) = Bessel's function. For the present it is 

m c 

sufficient to suppose < | e \ < M = finite constant. Considering fixed, for the 
moment, the upper limit of the outer integral, we can write 

$( 6 ) = c J°° n^T de f e^ ds. (IT) 

It is easy to see that 

(I) = real part of [^~] , (HI) 

because from (II') we get 

de J L z J 

the real part of which, after putting z — cs and remembering s= ~-, gives (I). 

* Sc. Pap., Vol. IV, p. 260. 
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We are going to transform the integral (II), and for this reason we take 
the well-known formula* existing for BessePs functions, 



J m (as)'J.(b8)= ^^ r^.(«Va' + 6'-2a6eoB»)^, d 
2"V*n(n— i) J o (Va 2 + 6 2 - 2 ab cos <p) n 

Let us put here n = l, a = b = c=-jj- • further remembering that Vl — cos<|> 
= V2sin^, then putting % = 4s and sin 4' = ?/, we have 

L « J V^t-n(i) Jo « 9 v a s 

Multiplying this equation by e i6 * and integrating from to oo with respect 
to s, we have 

Putting in the integral on the left cs = z and remembering that - = 2 e, then 
integrating on both sides from to e, we get 

»«=*XTPf i ] v " ,is ' i « 



_2c 
Let us consider 



r-Ji(2 C y5) fltti ^ 

Ja S 



Putting here 2cys = t and considering y as constant during the integration, 
we get 

Butt 



7 1= f ^iA^ e «-¥-<^. 

Jo t 



C x il^l- e ixt -dt = e* arcsin * for 0<»<1. 



e/x(Q . 

So we will have 

7 1 =e iarcsin (|), (1) 

* Nielsen, "Zylinder Funktionen, p. 183. f Nielsen, loc. cit., p. 197. 
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if 0< — <1. * In this inequality lies the definition of e. Since y varies from 
to 1 , it is evident that the domain of variation for s is 

0<* = £<l.t 
So from (II") we will have, after substituting (1), 



or 



= vM(V) So dF S!v Vl-/[cos(arc sin ±) + i sin (arc smj)]dy. (II"') 

Now we have | 

sin (arc sin x) = ( — 1)" sin [(— l) n |arc sin x\] 

and 

cos (arc sin x) = ( — l) n cos [(—1)" |arc sin x\]. 

The first of these two functions does not change its sign when n is odd or even; 
the second changes. But we have 

sin ( | arc sin x \ ) = x 
and 

Therefore 
and 

So in our case 



and 



And (II'") becomes 



cos ( I arc sin x \ ) = Vl — # 2 - 
sin (arc sin x) = x 



cos 


(arc 


sin x) = 


-.±y/l- 


-x\ 




sin ( arc 

arc sin - 

V yi 


sin 

1 = 


± 


H 




COS | 


^- 


e 2 

y 




1 «* 



* See Genocchi, "Diff.-Int. Rechnung," p. 170, for lim -1 >.0. 

e=0 \VJ 

t Evidently, Jf = l. 

i Laska, "Formelsammlung," p. 2. 



■ 0<£<1. 
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We will deal separately with the following two integrals : 

l s = f'YT^y'dy, 
After rearrangement we have 

h = is f\l^y"Jl-^dy, 
which is an elliptical integral. Further, we can write 

r -iJ r 1 d y -d^. 1 \ C y%d y + 1 c x y id y~\ m 

where y/W = yj (1 — y 2 ) (l — ^j . Or briefly 

the meaning of A, B and G being clear. 

Now we can reduce the C by the well-known reduction formula * 



-2(n-2)b cy^*Ay 

J \/W 



i+i 



r 1 

Putting here w=4, a = l, b = ^ — , c= - 2 and taking into consideration the 

limits of our integral, we get 

and so (2) becomes 

h = u[lA-l(l + l) B ] (2<) 

But we can evaluate A and £. In 



yr^i_; 






* Legendre, " Traits des f onetions elliptiques," T. 1, p. 4. 



102 Steic: Note on the Theory of Optical Images. 

put y = ex and we get 

A _ C\ & X 

~ Jo VI — e 2 * 2 Vl — a 2 ' 
which is equal* to 

^ = e (if + *#'), 

where the K and IT' are complete elliptic integrals of the first species, with 
the moduli s and Vl — e 2 respectively. Further, 



B 






becomes by the transformation y = e x 

f i j£dx r r 1 r|] 

e J Vi-6 2 ^ 2 Vi-« 2 IJo^JlJ' 
But we can writef 

Z? = e [*#' + # — #], 

where i? and 2?' are complete elliptic integrals of the second species, with 
the moduli e and Vl — £ 2 respectively. 
So (2) can be written 

h = ie]^e{K + iK')-±(l + ~^- e {iE' + K-E)^. 

We get immediately 

I 3 =f\T=y~*dy = ?. 

Substituting these values of I 2 and l z in (II IV ), we obtain 

vv ; V«n(i) Jo L 13 

Separating the real and imaginary parts, 

* Durege, "Theorie d. ellipt. Funktionen," p. 91. f Glaisher, Qiittrt. Journ., Vol. XX, p. 315. 
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Now remembering (II) and (III) and taking into consideration II (i) = 
^, we get 

1= C Ri( c *) T cos bsds= =r^- [ (1 + e 2 ) E'-2 e*K'] 

and /' = pr-M^fl] 2 s inbsds=ce±i?-[(l+e i )E-(l-i?)K], 
Jo L s J oft 

where 0<(e = ^— )<1. 

In (a) we cannot retain the double sign, and we are going to decide which 
sign is to be used. As for the integration of <p (e) , we observe that we can use 
the following formulas: 



(a) 



(»+2) fe n Ede = s n+1 E + fe n Kde 
and (n + 2) 2 }e n+2 Kds = s n+1 E — (n+2)e n+1 (l — e 2 )K+ (n + l)*fe n Kde.} 

If in (a) we put e = 0, or, what is the same, & = 0, we get 

i im rr^T cos6s , s= c\^r\ds=^, 

i=0 Jo L S J * ; o L 5 J 3 71 

since 



(3) 



t 

€=0 e=0 



n(i) = r(|)=-Y ? , lim#' = l, lime 2 Z' = 0. 



If we put r = jtt = l and k = cs, in the known formula* 

r-JAk^jjjjci ,, r(j)-r(y+y) . 

Jo ""' r+" aK 2"+"-r(r+i)r(^+i)-r( v +^+i) ' 



we have 



r Ri(M) T rf ,-+4c. 

■/o L s J 3 71 

— 1V ; by taking the limit for v = 0. 

- zys 2 — v 2 

We see that lim 7>0 ; on the other hand, we can show by simple expansion 

e=0 

that for all values of e the factor [(l + e 2 )E' — 2£K'\ in the expression of I 
is always positive. Therefore, I is positive for all values of e, and we must 
drop the negative sign from its expression. We will decide later on the choice 
of sign in /'. 

* Nielsen, loo. oit., p. 194. f Gray-Matthews, " Bessel's Functions," p. 201. 
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It remains to investigate the case for e = ^- > 1. To distinguish this case 
from the previous one, we will let e = e', and we can write as before 

4>( f ') = c j* £ [hT^T ei2e " dz ds '> ( IV ) 

where it is sufficient to suppose |e'| > 1. Again, 

(I)=realpartof p^f) ]. (V) 

Omitting superfluous calculations, we can write immediately 

4,(6') = !? f ds' CyVT=?dy C Jl{2cy ' x) e ibx -dx. (IV) 



Now 



Jo x Jot 



if 2 c y - x = t and pr— = e'. But* 

C i±ifL e i * t dt = i(x-Vx T =l) 

Jot 

if x > 1. So we will have 



MHG)- 1 ] 



2/ ^y. 
if — >1. Since « varies between and 1, in order that e' should satisfy this 

y 

inequality for all values of y, its domain of variation generally must be 

e'>l.t 
So (IV) can be written 



♦ <0=*£f *#VX=5<(f -J£ -l)*. <iv-) 

We have to evaluate the integrals 
h = fy^fdy = -J 



and 



Is = f o y-Vl-y 2 'yj~-l'dy = e'f o Vl~y 2 -^l-^-dy. 

* Nielsen, loo. oit., p. 197. 

f— >.l could be satisfied, for |2/|< 1, by |e'|< 1; since our inequality must be satisfied for all val- 
ues of y, even for y = \, therefore we must take |e'|,2. 1. 
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Introducing — 4 = e" 2 <l and applying the previous transformations and reduc- 
tions, we get finally 

where the modulus of K and E is— T = e" < 1. 

e 

Substituting I 2 and 7 3 in (IV" ) , we get 

,, ,. Ac re' Sn , e' 2 — 1 ^ Z+lpl,, 

We are able to evaluate <p(e') with the aid of formulae (3), but we are not in- 
terested in the explicit value of fy(e'); we are interested in ^, , ' , in the 
equality 

^e' 7t L 4 3 3 J 

The real part of this expression is equal to zero; therefore, by (V), 

T = jTi^lT cos bsds ==. 

and J'= fT^^Tsin&sds = ce'+^[(e' 2 -l)li:-(<?' 2 +l)#], 

»/0 L s J olt 

where (V = -^— J > 1 and the modulus in Z" and JB7 is — . 

If, in Z', we put e'= 1, we get 

limZ'=c 8c 



(P) 



e'=l o7l 

On the other hand we have for V in (a) for e = 1 

hm 7' = c ± -=— . 

e=i o7t 

Since the integrand in I' is a continuous function, we must have 

lim r'=lim/'; 

6'=1 6=1 

further, since, in /', (1 + e*)E — (1 — e 2 ) K > for all values of s, we conclude 
that we must drop the positive sign in the expression of T. 



(a') 
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So we can finally write the formulae (a), 
1 = JT^lMTcos bsds = ~ [ (1 + e*)E'- 2e 2 K'\, 

I' = J^pMMT gin fcsds = ce--^- [ (1 + e*)E + (1 - e 2 )#], 

where 0<(£=-^— J<1. 

It is interesting to investigate whether some of the functions I, I', Y, Y' 
show maxima or minima. We take T and write, after putting c = 1, the con- 
dition for maxima and minima 

^ = 1 - 4~ [~(1 + f 2 ) ~ + 2e# - (1 - 6 2 ) df + 2e^l = 0. 
as 3n [_ «f «f J 

Or, since* 
and 



we have 



Hence 



"=-?<*-*>■ 






!B =^ = 0.7854 = |«{l-(^)V + 3( 2 -ij)V-5(^)V + ....}. 

By trial we find that e = 0.543 approximately satisfies this equation. That 
this is a maximum point, we see immediately from 

f^= — — [2#(0.543) — iT(0.543)]. 
as 7t 

From tablesf we pick out the values of E (0.543) and If (0.543) and get 

~ =— — [2X1.445 — 1.712] <0. 
at % 

Similarly we find that neither I nor Y' has any maxima or minima. 

What is the value of litn J'? We put c = 1 and — r = v\ in J', and get 

6' =00 6 

iimr'=iim- 3 T>? 2 -4-i(H->? 2 )^('7)-(i-'7 2 )^(>7)sl=iim^-. 

* Laska, Formelsammhmg, p. 329. f Jahnke, " Funkt. Tabl./' p. 68. 
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If we substitute yi = 0, we get the indeterminate form j- . By repeated appli- 
cation of the rules for the determination of such forms, we get finally 

lim r'=lim Z'=0, 

ri=0 e'=oo 

which shows that the curve representing the function Y' approaches asymp- 
totically the e'-axis. 

The following tables contain numerical values for different values of 
-e and e' : 



e = sin0 


*<«) 


''(e) 


—: = sin 

6 


r'(«') 


0=0° 


4c 
3n Xl 





= 0° 





5° 


" X 0.960 


0.079c 


5° 


0.011c 


10° 


" X 0.882 


. 140c 


10° 


0.012c 


15° 


" X 0.776 


0.192c 


15° 


0.023c 


20° 


" X 0.667 


0.228c 


20° 


0.038c 


25° 


" X 0.545 


0.249c 


25° 


0.047c 


30° 


" X 0.435 


0.258c- 


30° 


0.060c 


40° 


" X 0.247 


0.252c 


40° 


0.089c 


50° 


" X 0.112 


0.225c 


50° 


0.106c 


60° 


" X 0.040 


0.194c 


60° 


0.115c 


70° 


" X 0.007 


0.169c 


70° 


0.130c 


80° 


" X 0.001 


0.154c 


80° 


0.146c 


90° 





0.150c 


90° 


0.150c 



Graphically : 




In the figure we take, for unity of the ordinates, ^~ = 60 mm. in the case 

OTt 

of l(t) ; 200 mm. in the case of I'(s) and T'(e'). The line of abscissas con- 
taining e = sin $ refers to 7(e) and I'(s) ; the one containing e'= sin 6 refers 
to Y'(e'). 
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We now come to the end of our investigation. In our problem we are 
concerned only with the integral 



for 



or explicitly 






-A- = *L [ (1 + S *)E'- 2e*K']. 



To every value of b will correspond a value of — — . We get the greatest value 

A b 

for — when b = : — - = -= — . The smallest value of — is for b = 2c : -— = CL 

A b A OTt A b A 2c 

In this ease we have no image. For details see Mandelstam's work. 



